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Abstract. In the context of communication networks, the framework of stochastic event graphs 
allows a modeling of control mechanisms induced by the communication protocol and an analysis 
of its performances. We concentrate on the logarithmic tail asymptotics of the stationary 
response time for a class of networks that admit a representation as (max,plus)-linear systems 
in a random medium. We are able to derive analytic results when the distribution of the holding 
times are light-tailed. We show that the lack of independence may lead in dimension bigger 
than one to non-trivial effects in the asymptotics of the sojourn time. We also study in detail 
a simple queueing network with multipath routing. 

1. Introduction 

In this paper, we study tail asymptotics of the form 

lim -logP(Z > x) = -6*, 

x^oo x 

where the random variable Z corresponds to a "global" state variable associated to a (max,plus)- 
linear system. We only deal with light-tailed distributions, i.e. distribution functions that are 
decaying exponentially fast. The simplest example of random variable covered by our results 
is the stationary waiting time in a FIFO GI/GI/1 queue. This case has been extensively 
studied in the literature and much finer estimates are available, see the complementary works 
( jlglehart 1972] ) and (|Pakes 1975ft . 



In recent years, there has been some interest in extending this result to networks of queues. 



( Chang 1995 ) considers an intree network and uses large deviations techniques to show that the 
queue length distributions have an exponentially decaying tail. (Gane sh and Anantharam 1996) 
obtains the decay rate of the tail distribution for two exponential server queues in series fed by re- 



newal arrivals. In (Bertsimas, Paschalidis, and Tsitsiklis 1998), the decay rate of the stationary 



waiting time and queue length distributions at each node is computed in an acyclic network in 
the context of quite general arrival and service processes. Literature on large deviations of queue- 
ing networks with feedback is rare and confined to the setting of networks described by finite- 



dimensional Markov processes, see (Dupuis and Ellis 1995), (Dupuis, Ellis, and Weiss 1991) and 



the recent works (Ignatiouk-Robert 2000), ( Ignatiouk- Robert 2005 ) . Moreover, these works con- 



centrate on local large deviations and cannot handle the large deviations of the network in its 
stationary regime. The large deviations asymptotics of queueing systems are difficult to analyze 
because they are dynamical systems with discontinuities. To the best of our knowledge, there 
is no rigorous result on the large deviations of non-exponential networks with feedback in their 
stationary regime. 

In this paper, we consider a class of networks that admit a (max,plus)-linear representation. 
This class contains the stochastic event graphs (which can be used to model window-based 
congestion control mechanism like TCP) and hence our results give the tail asymptotics of the 
steady state end-to-end response times of these networks. We should stress that the results 
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of this paper are not restricted to this sub-class and we give an example of a network with 
multipath routing that is covered by our framework. 

From a mathematical point of view, we study Z the stationary solution of a (max,plus)-linear 
recursion. Precise results concerning large deviations of products of random topical operators 
have been obtained in (Toomey 2002). However very restrictive conditions are required on the 



coefficients of the matrix and only the irreducible case is studied in ( Toomey 2002| ). Here we 



do not assume these requirements to be fulfilled and we show that under mild assumptions 
on the matrix structure, the tail behavior of Z is explicitly given and can be computed (or 
approximated) in practical cases. 

In the next section, we first give the general (max, plus) framework, with some examples of 
queueing networks. Then we give the stochastic assumptions and the tail asymptotics of the 
stationary solution of the (max,plus)-linear recursion is derived in Theorem [T] which clearly 
extends the case of the single server queue. Theorem [2] gives a more explicit form of the rate of 
exponential decay. 

In Section El we study two queueing applications. First we consider a system of two queues 
in tandem and show that when the service times at both queues are identical, then depending 
on the value of the intensity of the arrival process there is a phase transition in the behavior of 
the network reaching a large end-to-end delay (Proposition [2]). Then we study in detail a simple 
example of queueing networks with resequencing. Multipath routing has recently received some 
attention in the context of both wired and wireless communication networks. By sending data 
packets along different paths, multipath routing can potentially help balance the traffic load and 
reduce congestion levels in the network, thereby resulting in lower end-to-end delay. We show 
how our framework can model such mechanisms and give analytical insights. 

Sections (H [5] contain the proofs of Theorem Q] and Theorem [2] respectively. We give some 
further directions of research in the conclusion. 



2. Tails for discrete event systems 



In this paper we consider open systems with a single input marked point process N = 
{(T n , An, -B n )}_ oo<n<0O , where in a queueing context the sequence {T n } describes the arrival 
times of customer in the network and {A n ,B n } carries the information related to the n-th 
customer (like its service time at the different stations, routing decisions). We give a precise 
description of the dynamic of the system in the next section and of the random variable Z for 
which we derive the logarithmic tail asymptotics. 



2.1. (Max, plus)-linear systems. The (max, plus) semi-ring M max is the set 1R U { — oo}, 
equipped with max, written additively (i.e., a © b = max(a, 6)) and the usual sum, written 
multiplicatively (i.e., a<8>6 = a + 6). The zero element is — oo. For matrices of appropriate sizes, 
we define (A®B)W> = A^®B^ := max(A^\ B^)) an d (A®B)M> = ($ k A^®B^ := 
maxfc (A^'V + B^). By convention if A is a matrix and c <E K max , then (A®c)^ ■= AM®c 

Let s be an arbitrary fixed natural number. We assume that we are given with a sequence of 
matrices with non-negative coefficients: A n of size s x s and B n of size s x 1. To the sequences 
{A n } n , {B n } n , and {T n } n , we associate the following (max, plus)-linear recurrence: 

(1) X n +i = A n+ i <g> X n © B n+ x (8)T n+ i, 
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where {X n , n £ Z} is a sequence of state variables of dimension s. In Examples [TJ [21 we derive 
the explicit form of this recursion for the single server queue, queues in tandem and a fork join 
system. We refer to these examples to get an interpretation of the various quantities. 

The stationary solution of this equation is constructed as follows. We write 

(2) Y[m,n] ■= D [k+l,n] ® B k <g> T k , 

m<k<n 

where for k < n, D\k+i,ri\ = ®j=n A? = A n ® ■ ■ ■ ® A k+ i and D[ ra+l n ] = E, the identity matrix 
(the matrix with all its diagonal elements equal to and all its non-diagonal elements equal to 
— oo). It is readily checked that Y[ m m i = B m <g> T m , and for all n > m, 

Y[m,n+1] = ® Y [m,n] © B n+ i ® T n+ \. 

In view of ^j, the sequence {X[_ nj o]} is non-decreasing in n, so that we can define the sta- 
tionary solution of dU), 

*?-oo,o] := lim Y [-n,o] ^ 00 • 



We define the stationary maximal dater by 

(-oo,0] 



(3) < Z := = r<? - To < oo. 



l<i<s 

The following writing for the stationary maximal dater shows the similitude with the traditional 
stationary workload of a single server queue: 

(4) Z = sup{S n -(T -T n )}, 

n<0 

where the process {S n } n <o is defined by 

(5) S n := (D [n+lfi] ®B n ) {i) . 

l<i<s 

Example 1. Consider a FIFO single server queue where T n is the arrival time of the n-th 
customer and a n is its service time. Equation ([1]) is then the standard Lindley's recursion, 

Xn+i = max(«Y n + (T n+ i,T n+ x + <J n +i) 
= a n+ i <g> X n o- n+ i (8) T n+1 . 

The interpretation of <Y n is the departure time of the n-th customer from the queue. Note that 
in this case, we have S n = Y2i=n °* an< ^ % is the stationary workload. 



Example 2. Consider now a system of two queues in tandem, where T n is the arrival time 
of the n-th customer in 
Equation ([T]) is given by 



(i) 

of the n-th customer in the system and cr n is its service time at queue i, for i = 1,2. Then 



X n+l 1 












(1(32) 
\ "n+1 




n+lj 



where we used the shorthand notation = ^+i ® c n +i- In this case <f„ is the departure 

(2) 

time of the n-th customer from the first queue and X n from the second queue. Hence Z is the 
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stationary end-to-end delay of the network. Note that in this case, we have 

£ 



(6) 



Sn 



— — i=n j=l 



Example 3. Consider the standard fork and join system as depicted (with Petri net formalism) 

in Figure [TJ In this model, each time a packet (say k) finishes its service in node 1, there is 
one packet sent up and one packet sent down simultaneously. The 'up'-packet ('down '-packet) is 
then also the k-th packet for node 2 (for node 3 respectively). The k-th packet joins the queue of 
node 4 once both packets have left node 2 and 3 respectively. Each node is a standard ■ /G/l/oo 
queue. 




Figure 1. Fork and join model. 



(i) 

Let Xn denotes the departure time of the n-th packet from node %. We have the following 
equations: 

{T n+l ®X^)®a^ +1 , 



X. 



n+l 

(2) 
n+l 

(3) 
n+l 

(4) 



(e*4 2) )^ (2) 



n+l' 

(2) 
n+l' 



Ln+1 - rf>4 3 'iffi4 4) )®o. 

This system is linear in the (max, plus) semi-ring R maX ) and we can write the recursion ([1]) with 
the following matrices: 

\ 







— OO 


— OO 


— OO 




(1®2) 
On 


a {2) 


— OO 


— OO 




(1®3) 
On 


— OO 


a (3) 

On 


— OO 


V 


(10203) 
o"n 


On 


C7 (3) 
On 






A. 



where we used the shorthand notations, an 



t 



B„ 



0". 



(0 



V o" n 



(1(8,2) 
On 

(1®3) 
On 

(l®2e3) 



(0 



In this case Z is the stationary end-to-end delay of the network. 



We refer to ( |Baccelli, Cohen, Olsder, and Quadrat 1992 ) for other examples of (max,plus)- 
linear networks (see also (Lelarge 2005) for an example showing how to model window control 
mechanism) . 
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2.2. Tail asymptotics for the stationary solution. We first need to give some condi- 
tions that ensure the stability of the system, i.e. that the limit ([3]) is finite. Hence we 
introduce the following assumption in order to apply first-order Theorems of Section 7.3 of 
dBaccelli, Cohen, Olsder, and Quadrat 1992[ ): 

Assumption (ST). (Structure of (A n , B n )) 

The random variables A n ,B n have a fixed structure, i.e. each entry of A n or B n is either a.s. 
equal to — oo or non-negative for all n. And each diagonal entry of A n is non-negative. 

Assumption (IA). (independence assumption) 

We suppose that the sequences {(A n , B n )} n and {r n := T n+ \ — T n } n are mutually independent 
and each of them consists of i.i.d. random variables with finite means. 

This assumption implies a law of large number for {S- n } defined in ([5]), namely, 

(7) — - — > n _>oo 7 both a.s. and in L%, 

n 

where 7 is a constant referred to as the top Lyapunov exponent of the sequence {A n } see 
Theorems 7.27 and 7.36 in dBaccelli, Cohen, Olsder, and Quadrat 1992[ ). 

Assumption (S). (stability) 
We assume that 7 < E[n] = a. 

We have that under [(I A ) | and | ( S ) | the maximal dater Z defined in ([3]) is almost surely finite. 

We denote by the vector with all its entries equal to 0. 
Assumption (SP). (Separability) 
We assume that we have for all n, 

4®o = B n eo. 



This assumption ensures that for a solution X n of ([!]) with any initial condition: if X n < 
(g> T n+ \ then the process X n +i-, Xn+2, ■ ■ ■ does not depend on the past X n , X n -i, ■ ■ ■ Note that 
this assumption is clearly satisfied in the examples described above since we have A n ® = B n . 



We refer to Section 2.2.4 of (Lelarge 2005D for an example of network with A n ®0 = £? n ©0 7^ B n . 



In fact, Propositions 3 and 4 of ( |Lelarge 2005D show that any FIFO event graph with a single 
input fits into our framework. However this condition allows also to deal with some type of 
networks with (random) routing as described in Section 13.21 This property of separability can 
be made precise in a larger framework than (max,plus)-linear networks: this is the class of 
monotone separable networks introduced in (B accelli and Foss 19 95). 

We now give the stochastic assumptions that ensure that the random variable Z is light tailed: 

Assumption (LT). (Light-tailed) 

Let 



T] = sup h > 0, E p 1 '"' 



(i) 



< OO 



We assume that r/ > 0. 
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We will always assume that Assumptions (ST), (IA), (S), (SP), (LT) hold. We are now 
in position to state our main result. 

Theorem 1. The following limit exists as an extended real number: 



A s (0) 



1 



lim — log E 

n— >oo n 



(8) 
We have 

(9) 9* = sup{9 > 0, A s (0) + A T (-0) < 0} > 0, 

where At {9) = logE [e 9ri ] and the tail asymptotics of Z is given by, 

lim -logP(Z >x) = -9*. 

X^QO X 

In the case of the single server queue, we have clearly Ag(9) = logE[exp(#<7i)] and Theorem 
[1] is a standard result of queueing theory that goes back to the work (jCramer 19 38 ) and in a 
queuemg context to fllglehart 1972| ). We will give more comments on this theorem in Section [3l 



2.3. More detailed results. In this section we give a more explicit form for 9*. Without loss 
of generality, we may assume that the matrices A n have the following block structure: 

( A n (l, 1) | — oo | — oo | — oo \ 



A n (2,l) | A n (2,2) 



-oo 



-oo 



\ A n (d,l) | A n (d,2) | | A n (d,d) J 

where each A n (£,£) is an irreducible matrix. 

Theorem 2. Associated to the irreducible matrices {A n (£, £)}, we define the following function: 



Ai(0) 



lim — log E 

n— >oo n 



where the limit exists in 1U {oo} and is independent of u,v. Then we have 9* = minjry, 9 e } 
where the 9 E 's are defined as follows 

9 l = sup{# > 0, A £ (9) + A T {-9) < 0}. 

In the case of a single server queue with exponentially distributed service times, we have 
■q > 9* and this property remains valid for a large class of distributions. However, we show in 
the next section that as soon as we consider a network (i.e. with at least 2 nodes) then the 
parameter r] can play a role even with exponentially distributed service times. 

We first give a framework where r\ cannot play any role. Given a vector v = (v^\ . . . , v^ K '), 
we call a (max, plus) expression V a polynomial in v of unit maximum degree if it has the form 

where Kj C [1, K]. 

Corollary 1. If there exists a sequence of random variables {a n = (o4 , . . . , On )}n such that 
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(1) the components of a n are independent of each other; 

(2) for all i, there exists k such that A^'^ = o$ ; 

(3) each entry of A n that is not or — oo is a polynomial (in M max ) in a n of unit maximal 
degree. 

Then we have 9* = min{#^}. 

In a queueing context, the sequence of matrices {A n (£, £)} corresponds to a specific "com- 
ponent" of the network. It is well-known that the stability of such a network is constraint by 
the "slowest" component. Here we see that in a large deviations regime, if each component is 
independent of each other, then the "bad" behavior of the network is due to a "bottleneck" 
component (which is not necessarily the same as the "slowest" component in average). 

Remark 1. In the framework of last Corollary, the tail asymptotics for Z under heavy-tailed 
(more precisely subexponential) assumptions (i.e. when Assumption (LT) is not satisfied) has 
been derived in ( |Baccelli, Foss, and Lelarge 2004 ). In this case, the exact asymptotics (i.e. not in 



the logarithmic scale) are derived and the Lyapunov exponents of the sub-matrices {A n (£, £)} ap- 



pear. These exponents are known to be hard to compute ( |Blondel, Gaubert, and Tsitsiklis 2000[ ). 
Similarly in the light-tailed case, we see that the asymptotics in the logarithmic scale is given 
by the functions Aj>(9). The computation of these functions is not easy, in particular when the 
network has some feedback mechanism. More formally, we will see that this function is convex 
and its right-derivative at zero is exactly the Lyapunov exponent of the sub- matrices {A n (£, £)}, 
so that knowing the function Ai(9) allows to determine the Lyapunov exponent. In particular, 
at the level of generality considered in this paper, our result cannot be made in a more explicit 
form. We will see in the next section several examples for which the value of 9* has a simple 
expression in term of the parameters of the problem. It is interesting to note that, as opposed 
to the heavy-tailed case, exact tail asymptotics for Z under light-tailed assumptions seem to be 
out of reach in the general framework of (max,plus)-linear networks. 

Example 4. Going back to the fork and join system described in Example O we see that the 
irreducible matrices are one-dimensional and we have for £= 1,2,3, 



At(9) = lo g E 



Hence if 9 e denotes the exponential rate of decay for the single server queue fed by {T n , o4 }, 
then we have 9* = min{#^} in the case where each sequence of service times at each station are 
independent of each other. 



3. Queueing Applications 

3.1. The impact of dependence. In view of ([!]), Z is the supremum of a random process with 
negative drift and to make the connection with the existing literature, we state the following 
result (for a proof we refer to (Lelarge 2006) Corollary 3.2): 

Proposition 1. Under Assumptions (IA) and (S) and if 

(1) the sequence {S^ n /n} satisfies a large deviation principle (LDP) with a good rate function 

i; 

(2) there exists e > such that A$(9* + e) < oo ; 
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where 6* is defined as in (0|). Then we have 

(10) lim -logF(Z > x) = -6* = - inf 

x^co x a>0 a 

This kind of result has been extensively studied in the queueing literature and follows directly 



from the work (Duffy, Lewis, and Sullivan 2003). However, we see that considering the moment 



generating function instead of the rate function allows us to get a more general result than 
(fTUl) since we do not require the assumption on the tail (which is essential for (fTOj) to hold see 
(Duffy, Lewis, and Sullivan 2003)). Indeed this assumption ensures that the tail asymptotics 
of ¥(S n — (To — T n ) > nc) for a single n value cannot dominate those of P(M > x). In this 
case, equation (fTUj) has a nice interpretation: the natural drift of the process S n — (To — T n ) is 
(7 — a)n, where 7 — a < 0. The quantity 1(a) can be seen as the cost for changing the drift of 
this process to a > 0. Now in order to reach level x, this drift has to last for a time x/ot. Hence 
the total cost for reaching level x with drift a is xl(a) /a and the process naturally chooses the 
drift with the minimal associated cost. This can be made precise in some cases by a conditional 
limit theorem that characterizes the most likely path. 

In this section we show that this interpretation might be misleading in a queueing context. 
We consider a very simple example: a system of two queues in tandem. We assume that the 
sequence {(on , an^)} n is a sequence of i.i.d. random variables with max{E[oj ], E[cr^]} < a 
and E[exp#(<7^ + cj 2 ^)] finite in a neighborhood of the origin. 

If the service times of station 1 and station 2 are independent of each other, then the most 
likely cause of a given customer suffering a large delay is that a large number of its immediate 
predecessors require service times in excess of their inter-arrival times at one of the station. How- 
ever in the case where the service times are the same at both stations, we show that depending 
on the intensity of the arrival process A, two situations may occur: 

(1) if A < A c , then the most likely reason that a given customer suffers a large delay is that 
its own service time is large; 

(2) if A > A c , then the tail asymptotic of the end-to-end delay is the same as in the inde- 
pendent case. 

Let consider first the case where er„ and are independent. We are in the framework of 
Corollary [H Hence if we denote by 8^ the exponential rate of decay for the tail asymptotics 
of the station; 
then we have 



(i) 

of the stationary workload of a single server queue with arrival times T n and service times a„ 



lim -logP(Z >x) = -min(6» (1) ,6» (2) ). 

x-^oo x 

This result has been obtained in (Gancsh 1998). In words, we can say that the large deviation 
of the end-to-end delay in a system of two queues in tandem with independent service times is 
dominated by the "worst" one. 

Consider now the case where cr„ = for all n and the sequence {on } n is a sequence 
of i.i.d. random variables exponentially distributed with mean 1/fi. We assume also that the 
arrival process is Poisson with rate A < /i. Then a direct application of Theorem [2] gives, 

Proposition 2. In the previous framework, we have 

lim ilogP(Z > x) = -9*, 

x^oo x 
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with 

\>n/2 => e* = n-x. 

This proposition completes the result in (Gancsh 1998). For small values of A, the tail of the 
end-to-end delay is determined by the total service requirement of a single customer whereas 
when A > fJ,/2, it is the same as in the independent case. 

This shows that the behavior of tandems differs from that of a single server queue. In particu- 
lar (|Anantharam 1989]) shows that for GI/GI/1 queues, the build-up of large delays can happen 
in one of two ways. If the service times have exponential tails, then it involves a large number of 
customers (whose inter-arrival and service times differ from their mean values). This behavior 
is analogous to that of tandems where the service times are independent at each station or if the 
intensity of the arrival process is sufficiently large. If the service times do not have exponential 
tails, then large delays are caused by the arrival of a single customer with large service require- 
ment. In contrast, we see that a single customer can create large delays in tandems even under 
the assumption of exponential service times, if the intensity of the arrival is sufficiently low. 



3.2. A case of study: queueing network with resequencing. The aim of this section is 
to show that the results of this paper are not restricted to the class of event graphs and that 
our framework can deal with complex synchronization problems encountered in applications. 

In many distributed applications (e.g., remote computations, database manipulations, or data 
transmission over a computer network), information integrity requires that data exchanges be- 
tween different nodes of a system be performed in a specific order. However, due to random 
delays over different paths in a system, the packets or updates may arrive at the receiver in a 
different order than their chronological order. In such 3j CclS6 ? Si buffer (with infinite capacity) at 
the receiver has to store disordered packets temporarily. There is an extensive literature on rese- 
quencing problem and we refer the interested reader to the survey (Bac celli and Makowski 1 989) 
(see also the more recent work (|Jean-Marie and~G iin 1993j)). 

We consider a simple queueing model of disordering, namely a set of 2 parallel single server 
queueing station .|GJ|1 with renewal arrivals under probabilistic state-independent routing. 
This model constitutes an ersatz of the very complex situation one seek investigate. While the 
details of any protocol have been eliminated, the essence of network behavior (i.e. disordering) 
is preserved. In the sequel, we shall thus consider the model described in Figure [2 We assume 




node 4 

□ DP O ~ 



Figure 2. Resequencing problem. 



that packets arrive in the first queue according to a renewal process {T n }. We will model the 
routing at node 1 by a Bernoulli routing: with probability p (resp. (1-p)) the packet is sent up 
to node 2 (resp. down to node 3). Once packet k reaches the receiver, it leaves the system if all 
packets j with j < k have already left the system. Otherwise it stays in the resequencing buffer, 
where it waits for the packets with number less than k. 
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This model is similar to the standard fork and join system described in Example El Even if 
the routing mechanism is not the same here, we now show how we can adapt Example El In 
order to model the desired routing mechanism we will use the idea of clones, i.e., packets that 
behave like real packets except that they never require any service time: their service time is 
null. Suppose that the real route of packet k is up. Then at the end of its service in the first 

(3) 

node, a clone is sent to node 3. Since aY = 0, the departure time of packet k from node 3 is 



*f=rnaxO*f, 

In both cases the "real" packet k joins the queue of node 4 once "real" packet k 
it (and not before). In particular packets are ordered when they leave node 4. 



(3) 

Xu_i )• Similarly, if the real route of packet k is down, then a clone is sent up. 

1 has joined 



Let {Cn = (C 



n \ ■ ■ ■ Cn^)}n be a sequence of i.i.d. mutually independent random variables 
with finite mean and such that ®i ^ ex P(^Ci^ < 00 f° r 9 in a neighborhood of the origin. Let 
{r n }n£Z be a sequence of i.i.d. random variables, independent of everything else, with values in 
{2, 3}. We write P(r n = 2) = 1 — P(r n = 3) =: p, and assume that < p < 1. In order to apply 
our idea of clones, we consider the (max,plus)-recursion defined in Example El and we define 



0- 



(i) 



c 



(i) a m 



C (2) 1 



{r„=2}> 



a 



(3) 



C (3) 1 



{fn=3}) 



a, 



(4) 



0. 



We see that our system satisfies the assumptions of Theorem [2j In this case Z is the end- 
to-end sojourn time of a packet (taking into account the resequencing delay). The following 
proposition is a direct application of Theorem El In particular, note that r\ > 9 e , hence we have 



Proposition 3. In the previous framework, we have 

1 



lim -logP(Z > x) 



-9* < 0, 



where 6* = min{9 i } and the 9 i 's are defined as follows 

6 e = sup{9>0, Ai(8) + A T (- 

with 



<0}, 



Ai(0) 
A 2 (0) 
A 3 (0) 



logE 



in 



log [pK 



(2) 



log ((l-p)E 



+ 1 -p 



Ml 



(3) 



+ P ■ 



Then it is possible to make some optimizations. Assume to simplify that = for all n 

(2) (3) 

and that the sequence {Q, } and {Cn } are independent sequences of i.i.d. random variables 
exponentially distributed with respective mean 1/^2 and I//U3. We assume also that the arrival 
process is Poisson with rate A > max(^t2>^3); i- e - we cannot send all packets to one node. In 
this case, we have to load-balance the traffic in order to have a stable system. We have to chose 
p such that A < mm(fi2/p, ^3/(1 — ?)}■ in order to do so efficiently, we would like to maximize 
9*. Here, we have 82 = ^2 — Ap and #3 = fi^ — A(l — p), hence the optimal value of p is given by 



1 

P= 2 



A 



+ 1 and then, 9* 



l±2 + Ai3 - A 



In particular, in the symmetric case /xi = ^2, we find p = 1/2 which is in accordance with 
standard results of resequencing literature (1 Jean-Marie and Gun 1 993). 



TAIL ASYMPTOTICS FOR DISCRETE EVENT SYSTEMS 



11 



This model is certainly an oversimplification of the reality but it is intended to be also of 
pedagogical interest. We should stress that it could incorporate quite a few sophistication 
to enable to take into account general distributions, more than 2 paths or window control 
mechanism (where at any time, the number of packets in the network is upper bounded by the 
window size). In this cases, computations are much more complex and numerical approximations 
have to be made. 



4. Proof of Theorem Q] 

We first prove the existence of the moment generating function Ag given by the limit ([8]). 
Then we prove that 6* defined by ([9]) is positive and then we derive the tail asymptotics for Z. 



Note that under our assumptions we have for n > 0, 
(11) S n = d (D [lM ®B ) {i) , 

l<i<s 

where the equality is in distribution. Hence for the simplicity of notation, we will "inverse time" 
and index the processes by non-negative indexes. Hence S n is now given by the right-hand term 
of (jllj) and more generally we define for u < v, 



Ki<s 



4.1. Computation of the moment generating function. Note that thanks to Assumption 
(ST), we have with i* defined by S n = (-D[i, n ] ® -So) * > then 

Sn+l — ^-n+l ^ ^ n — '-"t 1. 

In particular the process {S n } is non-decreasing and we have the following rewriting: 

S luM =00 (%+M®^) W - 

l<i<s u<k<v 

The following lemma shows that the process has a subadditive property. 
Lemma 1. We have for n, m > 0, 

S n +m < S n + S'[ n+lin _|_ rn ] . 

Proof. We have by definition, 

D[k+i,n] ®B k <Q®S n , 

0<k<n 

hence by monotonicity, we have 

® D[k+i,n]®B k < A n +i (gO® S n , 

0<k<n 

D {k+1>n+1] ®B k < (B n+1 90)® S n , 

0<k<n 

iterating we get: 

D [k+1>n+m] ®B k < I D [j+hn+m] ® Bj J <g> S n . 

0<k<n \n+l<.j<n+m J 
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Hence we have 



D[k+l,n+m] 

0<fc<n+m 



Bk < 



D 



[j+l,n+m] 



S„ 



□ 



D 9S n 



,n+l<j<.n+m 

from which the lemma follows since Sj n +i )n + m i > 0. 
Lemma 2. The following limit 

A s (9) = lim ilogE 

n^oo n 

exists in MU {+00} /or a// 9 > 0. A^(.) is a proper convex function which is finite on the interval 
[0,r/). Moreover for all n and for all 9 <rj, we have E[exp 9S n ] < 00. 

Proof. In view of Lemma [U we have clearly the following subadditive property, for 9 > 0, 



logE 



< logE 



+ logE 



^9S„ 



and the existence of As follows and moreover, we have 



(12) 



A s (9) = inf ilogE 

n>l n 



The last part of the lemma follows from 

n 

i 

thanks to Assumption |(SP)| Then we have 



fc ' 



k=0 i 



n+l 



< E [e e5 " 



< 



(g)E 



hence we see that As(9) is finite for € [0, 7/) and infinite for 9 > n. Then the fact that Ag is a 
proper convex function follows from Lemma 2.3.9 of (|Dembo and Zeitouni 1998]) . □ 

4.2. Computation of 6*. 
Lemma 3. We have 9* > and 

A s (9) + A T (-9) < if 9 €(0,0*), 
A s (9) + A T (-9) >0 if 9>9*. 



Proof. In view of Assumption (S) and ([7]), we can choose n such that E[S* n ] < na. Then the 
function 9 1— ► logE[exp(0S' n )] is convex, continuous and differentiable on [0,7/). Hence we have 



1 



log E [e eSn ] + logE 



11 



,-6n 



E[S n 
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a) +o(6), 



which is less than zero for sufficiently small 9 > 0. By (|12p . we have for such 9, 



A s (9)+A T (-9) < - log E [e eSn ] + logE 



n 



-o-n 



< 0. 



Hence 9* is positive and the last part of the lemma follows from the convexity of As(9) + At(- 
□ 
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4.3. Tail asymptotics for Z. 
Lemma 4. We have, 

limsup-logP(Z > x) < -9*. 

x^oo X 

Proof. First note that by LemmaEl we have rj > 0* and for all < rj, we have E[exp(0SV t )] < 
oo. We denote = (Tq — T_ n ) for n > 0. For any < 9 < 9*, we have thanks to Chernoff's 
inequality, 



F(Z>x) = P(U n {S n -Sl>x}) 
< £>(S„ - S T n > x) 



n 

nA T (-6) 



e 



where the last series converge because we proved in Lemma [3] that for 9 < 9*, we have, 
1 



■logE[e e5 "l +logE 

n 



e -0n 



A 5 (0) +A T (-0) < 0. 

□ 

Lemma 5. We have, 

liminf ~logP(Z > x) > -9*. 

x^oo x 

Proof. We denote A(0) = A s (9) + A T (-9) and S T n = (T - T_ n ) as in previous proof. We 
only need to consider the case 9* < oo. We consider first the case where there exists 9 > 9* 
such that A(9) < oo. First note that the function 9 \— > A(9) is convex, hence the left-hand 
derivatives A'(9— ) and the right-hand derivatives A'(9+) exist for all 9 > 0. Moreover, we 
have A' (9— ) < A' (9+) and the function 9 i— > ^(A'(9— ) + A' (9+)) is non-decreasing, hence 
A'(9) = A'(9—) = A'(9+) except for 9 £ A, where A is at most countable. Since A(9) < oo for 
9 > 9*, we have A(0*) = and A'(0*+) > 0. To prove this, assume that A'(0*+) = 0. Take 
9 < 9*, thanks to Lemma| we have A(0) < 0. Choose e > such that < A(9* + e) < e|A(0)|. 
We have 

A(9* + e) -A(9) 
~e < 9* -0' 

which contradicts the convexity of A(0). Hence, we can find t < 9* + e such that 

< A(i), t i A. 
Note that these conditions imply t > 9* and A'(t) > A'(9*+) > 0. 

Thanks to Gartner-Ellis theorem (Theorem 2.3.6 in (|Dembo and Zeitouni 19981) ). we have 

(13) liminf - logP{S n - S T n > na)>- inf A*(x), 

where T is the set of exposed point of A* and A*(x) = sup 9>0 (9x — A(0)). Note that from the 
monotonicity of 9x — A(0) in x as 9 is fixed, we deduce that A* is non-decreasing. Moreover take 
a = A'(t), then A* (a) = ta — A(t) and a € T by Lemma 2.3.9 of (IDembo and Zeitouni 19 98). 

Given x > 0, define n = \x/a\. We have 

-logP(Z > x) > — log¥(S n - SZ > na), 
x na 
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taking the limit in x and n (while a = A'(i) is fixed) gives thanks to ([13 



1 , , to - A(t) 
liminf - logP(Z > x) > — > -t > 

x^oo x a 



€. 



We consider now the case where for all 9 > 9*, we have A(9) = oo, hence rj = 9*. Take K > 
and define S^ m] = S [n>m] \[f =n l{®jBf> < K) and Z K = sup n > (S , ^ n 0] - SI). We have clearly 
Z > Z K . It is easy to see that the proof of Lemma [2] is still valid (note that the subadditive 



property carries over to <S^ m j) and the following limit exists 



K K 



(0) 



lim — log E 

n— >oo n 



inf — log E 

n n 



OSr, 1 



Moreover thanks to the subadditive property of S (see Lemma [T]), we have S£ n ^ < Sff^ + 



+ 



5, 



A' 



jBi j) . Hence we have W(SZ n] < nil") = 1, so that Af (0) < Hence 



by the first part of the proof, we have 



lim inf — log ] 

x^oo x 



(Z K >x)> 



with 9 K = sup{6> > 0, Af (0) + A T (-9) < 0}. We now prove that 9 K -> rj as K tends to 
infinity which will conclude the proof. Note that for any fixed 9 > 0, the function Af (9) is 
nondecreasing in K and limx^oo Af (0) = A${9) < As(9). This directly implies that 9 K > 77. 
Take 9 > r], so that A 5 (0) = 00. If As(0) < 00, then for all K, we have Af (0) < As(0) < 00. 



But, we have Af (9) = inf n ±logE 



^tS'n 1 



so that there exists n such that 



E 



1, max(0 jj BJ- 



< e Af(0)+l < A<j(fl)+1 



but the left-hand side tends to infinity as K — > 00. Hence we proved that for all 9 > 77, we have 
Af (9) — > 00 as K — > 00. This implies that 9 K — > r/ as K — > 00. □ 



5. Proof of Theorem [2] and Corollary Q] 



We begin with a general result showing the existence of the function A^. Let {M n } be an 
i.i.d. sequence of irreducible aperiodic (max,plus)-matrices with fixed structure. We denote 



M, 



(hi) 



i [1)Ti] - (M n ®---^iivix i 
Lemma 6. For 9 > 0, the following limit exists in M U {+00} and is independent of i and j, 



A M (6) = lim ilogE 

n^oo n 



6M 



[l,n] 



Proof. 

We denote 



logE 



6M 



[l,n] 
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We have for 6 > 0, 



A 



(id), 
M I 



, n + m) 



logE 
logE 



0M, 



(id) 

[l,n+m] 



max e 



0M, 



> 



max <j log E 

k 



(i,k) 

[rc + l,n+m] g 
g [n+l,n + mj 



0M, 



(k,j) 

[l,n] 



max A 



M [ 



,m) + A 



M 



+ logE 



0M 



l,n] 



In particular for j = i, we have 

a£ j) ( 



n + m) > A 



M 1 



?,m)+A^ ) (0,n). 



Moreover thanks to the fixed structure assumption and the aperiodicity, there exists N such 



that for n > N, we have Mi^ > 



-oo for all i and j, hence Aw (9,n) > —oo and we have 



1 



lim -A&^n 



For arbitrary z and j, choose n,m > N and note that 



sup —K^\o,n) > — oo. 

n>iV " 



A 



M ^ 



?, n + ?rt) 



A^^n + m) 



> 
> 



A^(M) + A£%m), 
,n)+AS/ ) (0,m), 



where all terms are in R U {+oo}. Letting n — > oo while keeping m fixed, it follows that 



1 



lim -A ( if\e,n)= lim -A ( *f>(6,n). 



1 



n— >oo 77, 



n— >oo 77, 



□ 



Note that Assumption (ST) ensures that the matrices {A n (£, £)} are irreducible aperiodic with 
fixed structure. We now extend previous lemma to the sequence {A n } of reducible matrices. To 
do so, we first associate a graph Q = (V, £ ) to A n , as in Section 2.3 of ( |Baccelli, Cohen, Olsder, and Quadrat 1 992). 
Set V := {l,...,s}, which we abbreviate as [l,s]. An edge belongs to £ if and only if 

A»i > 0. Two nodes of V are said to belong to the same communication class if there is a 
directed path from the first to the second and another one from the second to the first. Let 
C\, . . . , Cd be the communication classes of Q and < the associated partial order, namely < C m 
if there is a path from any vertex in Ci to any vertex in C m . Without loss of generality, we assume 
that Ci <C m implies I < m; this is a notationally convenient restriction on the numbering of the 
communication classes. 

We use the following notation: 

• for any coordinate i € V, its communication class is denoted by [i], 

• for any coordinate i, the subset of coordinates j such that \j] < [i] is denoted by [< i]; 

• for any coordinate i, the subset of coordinates j such that [i] < [j] is denoted by [i <]; 

• for any coordinate i and j € [i <] , we write 

[i<j] ■= [i <]n[< j]. 

We now extend previous lemma. We introduce first some notations, 



A,(0) = A W (0) 



lim — log E 

n— >oo 77, 



6»D 



l,n] 
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which does not depend on i € Q as shown above. 
Lemma 7. For 9 € [0, rf), we have 



1 



lim — log E 

n— >oo 71 



(9D 



e l 1 .™] 

where the supremum over the empty set is — oo 



sup A[ fc ](0), 
fce[i<j] 



Proof. If [i < j] = 0, the result is obvious and if j € [i], the result follows from previous 
lemma. Hence we consider only the case: [i < j] D [i]. We denote 



A {i ' j) (9,n) = logE 



With the same argument as in previous lemma, we have for > 0, and for any [i] < Ce < [j], 
there exists k £ Cg and u, v > 1, 

A(* -J) (0, n + u + u) > A (i ' fc) (9,u) + A {k ' k) (9,n) + A^ (9,v), 

where each term is finite since 9 <rj. Hence we have by previous lemma: 

liminf SU p A m (0). 

n ke[i<j] 

(i j) 

For the upper bound, note that there exists u > 1 such that D^ J^ > — oo for all i < j. Consider 
first the case d = 2, i G C\ and j G C2, then we have 



A^' j \9,n) < logE 



9D^ a ] 9D<; a ' j) 



9D^' b) , eD\ h ' 3) 



maxe [Me +maxe [i.»-«]e [»-«+i,n] 

aeC 2 feeCi 



< log e A^(M e A«^(M + e .. v .,. ve . 



3 A( l ' a )(6»,«) p A( [l J)(6»,n) _j_ ^ p A( i . i, )(6»,n) p A( f, J)(6»,?i) 
v aGC 2 beCi 

Hence by Lemma 1.2.15 of (jDembo and Zeitouni 1998]) . we have 

A(*rf')(0,n) 



lim sup ■ 

We have clearly by induction that 



< max(Ai(0),A 2 (0)). 



A^)(0,n) 
lim sup - — < sup A[fci(0), 

n— >oo n ke[i<j] 



□ 



which concludes the proof. 

We now compute Ag(0) for a (max,plus)-linear system under the assumptions of Theorem [21 
Lemma 8. We have for 9 <G [0, rj) 

A s (9) =supA,(0). 



Proof. 

The lower bound follows directly from the following inequality: for all £, we have for i E Cg 



E e 



> E 



e [M 
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We now derive the upper bound. Note that A/. ® = © > B& , hence we have 

S n = © (D [k +l,n] ®B k ) {i) 



l<i<sO<k<n 
l<i,j<sO<k<n 



Hence we have 



E 



a SS n 



i,j k 



< n^E 



g [fe.n] 



e [0,n] 



and the lemma follows directly from Lemma 1.2.15 of (IDembo and Zeitouni 19981) . □ 

Theorem [2] follows directly form the fact that A<?(0) = oo as soon as 9 > r/ which follows from 
the lower bound Bfi' < S n . 

We now prove Corollary [TJ The following lemma implies that As (8) = sup^ Ai(9) for all 9 > 
and then Corollary [1] follows: 

Lemma 9. Under assumptions of Corollary {J\ there exists I € such that A^(0) = oo for 

all 6 > r/. 

Proof. We only need to consider the case r\ < oo. Take i such that sup{#, K[exp9B^] < 
oo} = r). By the condition on the entries of A n , we have 

K 



B^<(B 3 A^<Y,^ ] - 



k=l 



Hence there exists k such that K[exp(9a^)] = oo for 6 > ry and there exist j such that An 
Then t defined by Cg_ = [j] satisfies the property claimed in the lemma. 



(3,3) 



□ 



6. Conclusion 

We have shown that the distribution of the stationary solution of a (max, plus) recursion has 
an exponentially decaying tail and we gave an analytical way to compute the decay rate. 

We applied our results to different kinds of communication networks and exhibited quite non- 
standard behavior possible in high-dimension only. We also analyzed a queueing network with 
multipath routing and showed on a simple example how our analysis could help in the design of 
the routing decision depending on the characteristic of the traffic. 

We are currently working on some possible extensions of our work. Of particular interest 
would be a large deviations principle for the process {S n /n} introduced in ([5]). It would allow 
to give the most probable way for a large deviation of the maximal dater. Also, it should be 



possible to use the distributional Little's law QHaji and Newell 1971 [ ) to get asymptotics for the 
number of packets in the networks. 



In general, the characterization of the decay rate is given by the moment generating function 
which is not easy to compute, especially in the case of feedback. One practical question of 
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interest would be to find good ways to estimate this function from the statistics made on the 
traffic. 
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